finish
/clear
!
/filn,E353
/title,O-ring         !趨
rad=6.4               !Oλ뾶
th=1.28               !Oλں
/prep7                !ǰ
et,1,plane182         !ƽ浥Ԫ
et,2,CONTA172         !άӴԪ
et,3,TARGE169         !ӴĿ굥Ԫ
r,1                   !ʵ
mp,ex,1,2.148E+5      !ģ
mp,nuxy,1,0.3         !ɱ	      
mp,ex,2,210e6         !趨Ĳ2
mp,nuxy,2,0.3         !2Ĳɱ
TB,NLISO,1            !趨 
TBTEMP,0  
TBDATA,1,730,4000,300,800    

!Բģ
csys,1                !
k,1,rad-th            !ڱڹؼ
k,2,rad               !ڹؼ
kgen,36,1,2,1,,10     !ÿ10ȼؼ
*do,j,1,69,2          !趨ѭ
 a,j,j+1,j+3,j+2      !߳
*enddo                !ѭ
a,71,72,2,1           !ʣ߶ΣOλĴ
lesize,all,,,4        !߶Ϊ񻮷׼
allsel
amesh,all             !Oλ
nummrg,all            !±ŸԪر
numcmp,all            !ѹԪ

!·ģ
csys,0                !ѿϵ
rectng,-5*rad,5*rad,rad,5*rad     !Ϸ
rectng,-5*rad,5*rad,-rad,-5*rad   !·
!趨߶λֱ
lsel,s,line,,109,115,6
lesize,all,,,26,-0.2
lsel,s,line,,111,113,2
lesize,all,,,10
lsel,s,line,,110,116,6
lesize,all,,,6,2
lsel,s,line,,112,114,2
lesize,all,,,6,0.5
mat,2                  !2
asel,s,area,,37,38,1,1 !ѡ·
mshkey,1               !ӳ񻮷
amesh,all              !ֵԪ

!OλϷ֮Ӵ
lsel,s,line,,109       !ѡϷԲӴ߶
type,3                 !Ϊ3Ŀ굥Ԫ
real,1                 !ָʵ1
nsll,s,1               !ѡ߶ϵĽڵ
esln,s,0               !ѡڽڵϵĵԪ
esurf,all              !Ŀ굥Ԫ

allsel
csys,1                 !
lsel,s,loc,y,-225,-315 !ѡԲϷӴ߶
lsel,r,loc,x,rad
lsel,u,line,,109
type,2                 !Ϊ2ĽӴԪ
nsll,s,1               !ѡ߶ϵĽڵ
esln,s,0               !ѡڽڵϵĵԪ
esurf,all              !ӴԪ


!Oλ·֮Ӵ
et,4,CONTA172         !άӴԪ
et,5,TARGE169         !άĿ굥Ԫ
r,2                   !趨ʵ2
allsel
lsel,s,line,,115      !ѡ·ԲӴ߶
type,5                !Ϊ5Ŀ굥Ԫ
real,2                !ָʵ2
nsll,s,1              !ѡ߶ϵĽڵ
esln,s,0              !ѡڽڵϵĵԪ
esurf,all             !Ŀ굥Ԫ

allsel
csys,1                !ϵ
lsel,s,loc,y,225,315  !ѡԲ·Ӵ߶
lsel,r,loc,x,rad
lsel,u,line,,115
type,4                !Ϊ4ĽӴԪ
nsll,s,1              !ѡ߶ϵĽڵ
esln,s,0              !ѡڽڵϵĵԪ
esurf,all             !ӴԪ

!
/solu                 !
csys,1                !ϵ
nsel,s,loc,y,270      !ѡO·Ӵ
nsel,r,loc,x,rad
d,all,uy              !̶yλ

csys,0                !ѿϵ
nsel,s,loc,y,-5*rad   !ѡ·׶нڵ
d,all,all             !̶ȫλ
nsel,s,loc,y,5*rad    !ѡϷнڵ
d,all,ux              !̶xλ

outres,all,all        !趨ȫ
time,1                !趨ʱ䲽
nsubst,200            !趨Ӳ
nlgeom,on             !򿪴ο
lnsrch,on             !
nsel,s,loc,y,5*rad    !ѡϷȫڵ
d,all,uy,-0.24*rad    !ʩ12%ѹλ
allsel
solve                 !
*get,def1,node,166,u,y!ȡ166ڵyλƣΪdef1
nsel,s,loc,y,5*rad    !ѡϷȫڵ
d,all,uy,0            !λƻظ
allsel
solve                 !
*get,def2,node,166,u,y!ȡ166ڵyλ,Ϊdef2

disp1=(def2-def1)     !ûصdisp1
finish

/post1
*GET,ss,ACTIVE, ,SOLU,NCMSS        !ȡغɲ
*GET,ls,ACTIVE, ,SOLU,NCMLS        !ȡÿһغɲµӲ
*DIM,loads,ARRAY,ss,ls             !غ
*DIM,disps,ARRAY,ss,ls             !λ
*cfopen,spring,dat                 !springļ
*do,in,1,ls                        !ѭ,ȡÿһӲ½Ϣ
*do,jn,1,ss 
set,in,jn                          !ȡinغɲµjnӲϢ 
csys,1                             !ѡȡ·OӴ
lsel,s,loc,y,225,315
lsel,r,loc,x,rad
lsel,u,line,,115
nsll,s,1
FSUM,0,ALL                          !غ֮       
*GET,loads(jn,in),FSUM, , Item, Fy  !ȡyغֵloads
allsel
*get,disps(jn,in),node,166,u,y      !ȡ166Žڵλƣdisps
loads(jn,in)=loads(jn,in)           !鸳ֵ
disps(jn,in)=disps(jn,in)
*enddo  
*enddo 
!дļ
*do,kn,1,ls
*vwrite, loads(1,kn),disps(1,kn),
(1X,'    ',31F16.9,'    ',31F16.9)
*enddo  
*do,kn,2,ls
*vwrite, loads(2,kn),disps(2,kn),
(1X,'    ',31F16.9,'    ',31F16.9)
*enddo 